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Summary 

Resonance frequencies can provide useful information on tine deformation occurring during fractur- 
ing experiments or CO2 management, complementary to the microseismic event distribution. An 
accurate time-frequency representation is of crucial importance prior to interpreting the cause of res- 
onance frequencies during microseismic experiments. The popular methods of Short-Time Fourier 
Transform (STFT) and wavelet analysis have limitations in representing close frequencies and deal- 
ing with fast varying instantaneous frequencies and this is often the nature of microseismic signals. 
The synchrosqueezing transform (SST) is a promising tool to track these resonant frequencies and 
provide a detailed time-frequency representation. Here we apply the synchrosqueezing transform to 
microseismic signals and also show its potential to general seismic signal processing applications. 

Introduction 

Traditional time frequency representations, such as the Short-Time Fourier Transform (STFT) and 
the Wavelet Transform (WT) and special representations like Empirical Mode Decomposition (EMD), 
have limitations when signal components are not well separated in the time-frequency plane (Wu 



et al.|[20 11 ). Synchrosqueezing, first introduced in the context of speech signals JDaubechies and 



Maes[ ]T996) has shown to be an alternative to the EMD method paubechies et al.[[2011|^ , improving 



spectral resolution. 

On the other hand, these improvements in time-frequency representation can be applied to micro- 
seismic signals to increase the readability of the frequency spectrum ( [Auger and Flandrin| |1 995^ and 



identify independent components (| Thakur et~ar| |2012^ . An accurate representation of resonance 



frequencies could help gain a better understanding of the fracturing process, as it is challenging to 
determine what is exactly beneath the Earth's surface during the hydraulic fracturing process. 

In this paper we demonstrate the potentiality of the SST in time frequency representation of micro- 
seismic signals. First we describe the underlying modulation model for microseismic signals and then 
we validate the method with synthetic examples. Finally we apply SST to a real data sample with high 
amplitude noise. 

Theory 

Coherent signals in microseismic recordings can have both anthropogenic and natural sources. In 
both cases the signal has similar oscillatory characteristics and could be represented as the sum of 
individual time-varying harmonic components: 

K 

s(O = £a,(Ocos(0,(O) + t](O, (1) 

k=i 

where Au{t) is the instantaneous amplitude and fk{t) = ^j-,Qk is the instantaneous frequency (IF) of 
the resonator k. r]{t) represents the additive noise, including the contribution of environmental and 
acquisition sources. K stands for the maximum number of components in one signal. 
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Thus, microseismic signals show some similarities with the time-frequency representation (TFR) of 
speech signals (Daubechie s and Maes[ |1996). The number of harmonics or components in the 



signal (K) is random. They can appear at different time slot, with different amplitude {Ak{t)) and 
instantaneous frequency {/kit))- 

From CWT to SST 

The CWT of a signal s{t) is ( |Daubechiesl|1992[ ): 



W,{a,b) = -j=js{t)Y{ — )dt, (2) 

where v^* is the complex conjugate of the mother wavelet, b is the time shift applied to the mother 
wavelet which is also scaled by a. The CWT is simply the cross-correlation of the signal s{t) with a 
number of wavelets that are scaled and translated versions of the original mother wavelet. Ws{a,b) 
are the coefficients representing a concentrated time-frequency picture, which is used to extract the 
instantaneous frequencies (Daubechies et al., ,201 1 ) . 

The wavelet coefficients Ws{a,b), often spreads out around the scale dimension a, leading to a blurred 
projection in time-scale representation. Daubechie s and Maes| ( |1996| ) show that if the smear in the 



b-dimension can be neglected, then the instantaneous frequency C0s{a,b) can be computed as the 
derivative of the wavelet transform at any point {a,b) for which Ws{a,b) / 0: 

= db • 

The final step in the new time-frequency representation is to map the information from the time- 
scale plane to the time-frequency plane. Every point {b,a) is converted to {b,(Os{a,b)), this operation 
is called synchrosqueezing ( [Daubechies et al.[ |2011} . Since a and b are discrete values we can 



have a scaling step Aaj. = a^^i -a^ for any where Ws{a,b) is computed. Likewise, when mapping 
from the time-scale plane to the time-frequency plane {b,a) {b,Winst{a,b)), the synchrosqueezing 
transform Ts{w,b) is determined only at the centers co/ of the frequency range [w/ - Aa)/2,a)/ + A«/2], 
with Aw = ft)/ - coi i: 

U(Oub) = ^ £ W,(a,,%-3/2Aa,. (4) 



at : I ffl - 6); I < Aci)/2 



The above equation shows that the time-frequency representation of the signal s{t) is synchrosqueezed 
along the frequency (or scale) axis only ( |Li and Liang , 2012) . The synchrosqueezing transform reallo- 
cates the coefficients of the continuous wavelet transform to get a concentrated image over the time- 
frequency plane, from which the instantaneous frequencies are then extracted ( |Wu etal.[|2011f . This 



is an ultimate goal in microseismic signal analysis. The identified frequencies are used to describe 
their source mechanisms and eventually gain a better understanding of the reservoir deformation. 



Synthetic example 

Following the modulation model in Equation [T|we create a noiseless synthetic signal as the sum of 
the following components (the IF is in brackets): 

si{t) = O.3cos(107r0;[^^ = 5]; 

S2{t) = O.8cos(307r0;[/i^ = 15] 

si,{t) = 0.7cos(207r? + sin(7r?)); [/F = 10 + cos(7r?)/2]; 

-$4(0 = 0.4cos(667r; + sin(47r?));[/^' = 33 + 2cos(47rf)]; 
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The synthetic signal s{t) (Figure[T^) has two constant harmonics of 5 Hz {si{t)) and 15 Hz (^2(0) f''0'T^ 
time to 6 s, the 1 Hz component (5^3 {t))\s modulated by a 0.5 Hz sinusoid from 6 to 1 s and a 2nd 
modulated signal with a central frequency of 33 Hz (54(0) appears at time 4 s and vanishes at 7.8 s. 

Figure [ij)) shows the instantaneous frequencies for each component. The STFT is able to identify 
the four components but with a low resolution (see Figure [Tj:), especially when they overlap at 6 s. 
On the other hand, the SST (Figure is able to perfectly delineate each individual component and 
to resolve the instantaneous frequencies close to the theoretical value, including S3{t) and S4{t) which 
are amplitude modulated. 



Figure 1 Synthetic Example, a) Synthetic 
signal with four sinusoidal components, b) 
the instantaneous frequencies obtained as the 
derivative of each one of the independent com- 
ponents, c) The STFT of the synthetic sig- 
nal. A Manning window of 321 samples and 
50 % overlap was used. Note the smearing 
effect in the bandwidth. In d) the SST shows 
a sharper representation of the instantaneous 
frequencies. 



Application to real microseismic signals 

In this section we apply the SST to a real dataset. A two-stage microseismic experiment was recorded 
by twelve 3C geophones deployed in a vertical well. In both stages the sampling frequency was 4000 
Hz. The receivers and fluid injection are located approximately at the same depth. The SST will 
be challenged by typical problems encountered in microseismic experiments analysis, i.e, time and 
frequency resolutions, noise stability and capability to identify localized frequency variations. 

A segment of 5 minutes of data, presenting clear resonance frequencies as well as sharp and smooth 
changes in frequency content, was chosen to challenge the performance of the SST representation. 
The STFT will be used as a reference, see Figure |2] upper plot. The STFT was calculated using a 
window length of 2 s with 50 % overlap, involving a time resolution of the same order. Four main 
resonance frequencies are clearly visible at approximately 18, 31, 35 and 52 Hz (Figure |2]top). Two 
additional lines of lower amplitude are also visible at 60 Hz and between the lines at 31 and 35 Hz. 
At approximately 90 s, the amplitude of all lines becomes weaker, the one at 52 Hz almost vanishes. 
The line at 31 Hz abruptly drops down to 27 Hz while smooth changes are barely discernible on the 
18 Hz line. The fuzzy background is not only due to the high amplitude noise present in the data but 
also to spectral leakage introduced by the Fourier transform. 

Despite the high-amplitude noise, the resonance frequencies as well as their variations are all clearly 
visible on the SST representation (Figure |2] bottom). The SST is able to map smooth and sharp 
changes in the frequency lines. Compared with the STFT, the SST brings out the resonance fre- 
quencies more sharply, improving significantly the frequency resolution. The SST is then able to 
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distinguish between the different lines constituting the resonance frequency at 18 Hz while the STFT 
cannot. Last but not least, the SST determines the IFs at all times. Therefore, the time resolution of 
this method is not limited by the size of any window. 

For the SST we use a bump wavelet with a ratio central frequency to bandwidth of 50. The dis- 
cretization of the scales of CWT was 64. The SST produces a sharper time-frequency representation 
showing frequency components that were hidden in the STFT representation. 



150 
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Figure 2 Real example. Upper plot shows the 
STFT and the lower plot the SST The SST 
is able to delineate the spectral components 
some of them were missing in the STFT rep- 
resentation. 



Conclusions 



We have introduced the Synchrosqueezing transform with application to microseismic signal analysis. 
This new transform shows promising results in the identification of resonant components and hence 
in the explanation of microseismic phenomena. The enhancement in spectral resolution allows the 
separation of close spectral bands while the short-time Fourier transform smears out these frequency 
components. SST is therefore attractive for high-resolution time-frequency analysis of microseismic 
signals. 
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